#include <buola/mat.h>
#include <buola/mat/cgaussianprocess.h>

using namespace buola;

int main()
{
    mat::CVec_d lDataPoints{5,10,15};
    mat::CVec_d lDataValues{2,-2,3};
 
    mat::CGaussianProcess<double> lProcess(lDataPoints,lDataValues);

    for(double a=0.0;a<20.0;a+=0.1)
        msg_info() << lProcess.F(mat::CVec_d{a}) << "\n";
}